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14.  ABSTRACT 

The  effect  of  stress  ratio  on  the  statistical  aspects  of  small  fatigue  crack  growth  behavior  was  studied  in  a  duplex 
microstructure  of  Ti-6Al-2Sn-4Zr-6Mo  (Ti-6-2-4-6)  at  260°C  with  particular  emphasis  on  incorporating  small-crack  data 
into  probabilistic  life  prediction  and  the  influence  of  stress  ratio  on  probabilistic  lifetime  limits.  A  Focused  Ion  Beam 
(FIB)  was  used  to  machine  micro-notches  in  test  specimens,  which  served  as  crack-initiation  sites  and  enabled  the 
acquisition  of  multiple  small-crack  growth  data  sets  from  single  experiment.  Stress  ratios  of  -0.5,  0.05,  and  0.5  were 
employed,  and  small-crack  growth  was  monitored  by  the  acetate  replication  method.  Qualitatively,  change  in  the  stress 
ratio  produced  almost  negligible  influence  on  the  small-crack  growth  behavior  when  plotted  as  crack  growth  rate  (da/dN) 
versus  applied  stress  intensity  factor  range  (delta  K).  A  probabilistic  method  was  employed  to  represent  the  variability  in 
the  small-crack  growth  behavior  and  the  statistical  differences  with  respect  to  stress  ratio. 
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Abstract 

The  effect  of  stress  ratio  on  the  statistical  aspects  of  small  fatigue  crack  growth  behavior  was 
studied  in  a  duplex  microstructure  of  Ti-6Al-2Sn-4Zr-6Mo  (Ti-6-2-4-6)  at  260°C  with  particular 
emphasis  on  incorporating  small-crack  data  into  probabilistic  life  prediction  and  the  influence  of 
stress  ratio  on  probabilistic  lifetime  limits.  A  Focused  Ion  Beam  (FIB)  was  used  to  machine 
micro-notches  in  test  specimens,  which  served  as  crack-initiation  sites  and  enabled  the 
acquisition  of  multiple  small-crack  growth  data  sets  from  single  experiment.  Stress  ratios  of  -0.5, 
0.05,  and  0.5  were  employed,  and  small-crack  growth  was  monitored  by  the  acetate  replication 
method.  Qualitatively,  change  in  the  stress  ratio  produced  almost  negligible  influence  on  the 
small-crack  growth  behavior  when  plotted  as  crack  growth  rate  (da/dN)  versus  applied  stress 
intensity  factor  range  (AK).  A  probabilistic  method  was  employed  to  represent  the  variability  in 
the  small-crack  growth  behavior  and  the  statistical  differences  with  respect  to  stress  ratio.  In 
spite  of  qualitative  similarity,  differences  in  the  statistical  parameters  of  small-crack  growth  as  a 
function  of  stress  ratio  were  found  to  be  significant  enough  to  affect  life  prediction.  The  methods 
for  representation  and  probabilistic  treatment  of  small-crack  data  were  also  shown  to  be 
important  factors  in  incorporation  of  the  small-crack  regime  in  probabilistic  life  prediction. 

Keywords:  Small  fatigue  crack,  Stress  ratio,  Ti-6Al-2Sn-4Zr-6Mo,  Life  prediction 
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1. 


INTRODUCTION 


Accurate  characterization  and  understanding  of  the  effect  of  stress  ratio  (R)  on  fatigue 
crack  growth  are  important  from  several  life  prediction  considerations.  For  example,  under 
loading  profdes  of  practical  interest,  such  as  spectrum  loading,  the  R-effect  or  the  mean-stress 
effect  on  crack-growth  behavior  is  a  critical  input  in  lifetime  calculations  [1  -  3].  Similarly,  an 
accurate  representation  of  the  R  effect  is  required  in  problems  where  the  crack  grows  through 
regions  of  varying  stress  state,  such  as  a  weld  zone  [4,  5],  under  a  residual  stress  field  introduced 
by  surface  treatments  [6,  7],  and  crack  growth  from  geometric  features  [8]. 

The  effect  of  R  on  the  long-crack  growth  behavior  has  been  studied  extensively.  In  the 
long-crack  regime,  the  degree  of  the  R  effect  varies,  depending  on  the  material  [9,  10].  In 
metallic  materials,  the  effect  is  found  to  be  more  pronounced  in  the  near-threshold  crack-growth- 
rate  regime,  which  is  often  rationalized  in  terms  of  the  crack  closure  arguments  [11-15].  The 
crack  closure  effects  are  thought  to  decrease  with  increasing  R,  producing  a  shift  in  the  long- 
crack  growth  curve  towards  higher  growth  rates  [11  -  15]  when  plotted  in  terms  of  the  applied 
stress  intensity  factor  range,  AK.  It  has  been  shown  in  such  cases  that  the  R  effect  on  long-crack 
growth  largely  diminishes  when  data  are  plotted  in  terms  of  the  closure-corrected,  effective-AK 
(AKeff)  [11-15].  Although  crack  closure  appears  to  explain  the  R  effect  in  many  cases,  other 
mechanisms  have  also  been  shown  to  be  important  [16  -  20],  These  mechanisms  largely  stem 
from  the  recognition  that,  under  certain  loading  regimes,  AK  does  not  uniquely  represent  the 
driving  force  for  long-crack  growth  but  contributions  of  other  parameters  such  as  the  maximum 
stress  intensity  factor,  Kmax  need  to  be  represented  [16-20],  Thereby,  the  effect  of  R  on  long- 
crack  growth  has  also  been  explained  by  mechanisms  involving  static  fracture  modes  that  are 
promoted  by  Kmax  [16-20], 
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Small  fatigue  cracks,  which  by  one  definition  are  cracks  whose  dimensions  are  on  the 
order  of  the  microstructural  size  scale,  have  long  been  recognized  to  grow  below  the  long-crack 
growth  threshold  and  faster  than  a  long-crack  at  the  same  nominal  AK  [21-25].  As  such,  small- 
crack  growth  constitutes  a  vital  regime  in  fatigue  life  prediction.  This  is  not  only  because  a 
purely  long-crack  based  prediction  can  often  be  non-conservative,  especially  in  the  long-lifetime 
domain  [26,  27],  but  also  due  to  the  consideration  that  a  significant  fraction  of  lifetime  after 
crack  initiation  can  be  spent  in  the  small-crack  regime  [27,  28],  The  importance  of  appropriately 
accounting  for  the  small-crack  growth  regime  in  life  prediction  is  further  accentuated  by  the 
observation  that,  in  many  cases,  the  influence  of  a  variable  on  the  small-crack  behavior  is 
different  than  its  effect  on  the  long-crack  growth  [29,  30],  Additionally,  the  small-crack  behavior 
exhibits  a  stronger  influence  of  microstructure  [22,  24],  which  is  manifested  as  larger 
fluctuations  in  growth  rate  than  is  seen  in  the  long-crack  regime,  as  well  as  greater  variability  in 
growth  between  different  small-cracks  in  nominally  identical  material. 

The  above  factors  suggest  that  the  small-crack  regime  should  be  a  crucial  ingredient  in 
probabilistic  modeling  of  fatigue  lifetime,  and  recent  emphasis  on  physics-based  probabilistic 
life-prediction  methods  [31]  has  generated  a  renewed  focus  on  small-crack  growth.  However, 
although  the  small-crack  effects  are  considered  important  for  accurate  probabilistic  life 
prediction,  incorporation  of  this  regime  in  life  prediction  is  partly  hampered  by  the  very  nature  of 
small-crack  data,  which,  as  discussed  above,  is  marked  by  fluctuations  in  growth  rate  and 
variability  in  growth  between  different  cracks.  A  lack  of  sufficient  attention  in  the  literature  to 
quantitative  representation  of  the  small-crack  regime  and  the  effect  of  microstructure  and  loading 
variables,  such  as  R,  on  the  statistical  parameters  of  small  crack  growth  has  further  limited  an 
explicit  inclusion  of  small-crack  data  in  probabilistic  life  prediction  models. 
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From  a  purely  crack  closure  argument,  one  might  expect  an  almost  negligible  R  effect  on 
growth  rates  of  small  fatigue  cracks,  since  crack  closure  is  considered  largely  a  crack-wake 
effect,  whereas  the  wake  is  thought  to  not  have  fully  developed  in  the  small-crack  regime  [21, 
22],  A  lack  of  significant  influence  of  R  in  the  small-crack  regime  is  also  supported  by  a  recent 
study  on  Ti-6A1-4V  (Ti-6-4)  [32],  Other  researchers  have  shown  that  the  long-crack  behavior  at 
a  high  R,  where  the  crack  closure  effects  are  considered  essentially  negligible,  encloses  the 
small-crack  data  at  lower  R  values  [33,  34],  Life  prediction  models  for  fatigue  design,  based  on 
these  aspects  of  the  effect  of  R  on  small  vs.  long  crack  growth,  have  been  developed  [33,  34], 
Additionally,  models  of  plasticity-induced  crack-closure  that  correlate  the  effect  of  R  on  the 
long-crack  growth  and  account  for  the  small-crack  behavior  have  been  proposed  and  applied  to 
fatigue  life  prediction  [28].  These  approaches  have  provided  valuable  insights  into  the  limiting 
AK  conditions  and  methods  for  design  of  fracture  critical  structures.  However,  to  enable 
probabilistic  life  prediction  models,  also  needed  are  methods  that  address  the  quantitative 
representation  of  the  small-crack  growth  data  and  the  effect  of  R  (and  other  variables).  Some 
recent  studies  have  proposed  a  mechanism-based  probabilistic  model  of  life-limiting  fatigue 
behavior,  which  demonstrated  that  small-crack  growth  can  be  used  to  predict  the  life-limiting 
failure  distribution  [29,  35,  36].  From  a  probabilistic  life  prediction  perspective,  the  questions  of 
sensitivity  of  lifetime  to  variation  in  statistical  parameters  of  small-crack  growth  behavior  and 
the  data  reduction  method  are  considered  important,  and  this  paper  is  focused  on  these  aspects  of 
small-crack  growth.  In  this  regard,  the  effect  of  R  on  the  small-crack  growth  behavior,  both  in 
terms  of  the  average  growth  rates  as  well  as  the  variability  in  the  same,  in  the  a+p  titanium  alloy, 
Ti-6Al-2Sn-4Zr-6Mo  (Ti-6-2-4-6)  was  studied  at  260°C.  A  probabilistic  method  was  used  to 
represent  the  small-crack  growth  behavior,  which  was  applied  in  a  previously  proposed  life- 
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prediction  model  [29,  35,  36]  to  determine  the  effect  of  R  on  the  distribution  in  lifetime  and  the 
probabilistic  life  limits.  The  data  reduction  method  is  an  important  component  of  quantitative 
representation  of  the  small-crack  behavior,  where,  for  the  purpose  of  life  prediction,  the  effects 
of  the  measurement  error,  frequency  of  data  collection,  and  the  maximum  crack  length,  on  the 
modeled  da/dN  -  AK  curves  were  also  examined. 

2.  MATERIAL  AND  EXPERIMENTAL  PROCEDURE 

2.1  Material 

The  material  in  this  study  was  the  a+p  titanium  alloy,  Ti-6Al-2Sn-4Zr-6Mo.  The  alloy 
was  obtained  in  the  a+P  processed  and  heat  treated  condition,  which  produced  the  duplex 
microstructure  shown  in  Fig.  1.  The  microstructure  consisted  of  primary,  equiaxed  a  (ap) 
particles,  and  lamellar  secondary-a/p  colonies.  The  average  size  of  ap  was  about  4  pm  and  that 
of  colonies  was  about  15  pm.  The  volume  fractions  were  approximately  0.3  and  0.7, 
respectively,  of  ap  and  the  colony  phase.  Further  details  on  the  microstructure  are  provided  in 
[35].  The  tensile  behavior  of  this  microstructure  at  260  °C  is  shown  in  Fig.  2.  The  0.2%  yield 
strength  was  about  850  MPa  and  the  ultimate  tensile  strength  was  about  1010  MPa. 

2.2  Experimental  Procedure 

Small-crack  growth  testing :  The  small-crack  growth  tests  were  conducted  using  an  MTS 
servo-hydraulic  test  system  equipped  with  a  458  controller.  The  tests  were  performed  in  load 
control,  and  the  specimens  were  axially  loaded  at  a  frequency  of  20  Hz  at  260  °C  in  lab  air 
atmosphere.  Three  stress  ratios,  R  (omjn/amax)  =  -0.5,  0.05,  and  0.5,  were  used. 
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The  fatigue  specimens  had  a  round-bar  geometry  with  a  uniform  gauge  section  of  about  4 
mm  diameter  and  a  gauge  length  of  about  12.5  mm.  The  final  machining  step  was  low  stress 
grinding  (LSG).  Specimens  were  electropolished  to  remove  about  50  pm  of  the  surface  layer  in 
order  to  minimize  any  effect  of  LSG-induced  residual  stress  and  to  provide  a  uniformly  smooth 
surface  for  observation  and  measurement  of  small  cracks.  Nine  micro-notches  were  machined  in 
each  specimen  using  a  Focused  Ion  Beam  (FIB)  to  create  crack  initiation  sites  and  to  enable 
multiple  small-crack  measurements  from  each  specimen.  The  placement  scheme  of  FIB  micro¬ 
notches  was  similar  to  that  reported  in  [32],  The  micro-notches  were  arranged  in  3  columns,  each 
containing  3  notches,  where  the  columns  were  separated  by  120°  around  the  specimen 
circumference.  The  spacing  between  notches  in  a  given  column  was  3  mm.  Further,  the  three 
columns  were  staggered  by  1  mm  with  respect  to  each  other  so  that  only  one  notch  was  present 
on  any  given  plane.  An  example  of  a  FIB-machined  micro-notch  is  presented  in  Fig.  3.  In  the 
overview  fractograph  of  Fig.  3(a),  a  fatigue  crack  that  initiated  from  a  FIB  notch  in  one  of  the 
small-crack  growth  specimens  (tested  at  R  =  0.5)  is  shown  at  final  fracture.  This  figure  indicates 
that  the  crack  aspect  ratio  (i.e.,  depth/half  surface  arc  length)  approximately  equaled  1.  The 
higher  magnification  image  focused  at  the  origin  of  the  same  crack  is  presented  in  Fig.  3(b),  and 
shows  that  the  crack  emanated  from  the  FIB  notch.  The  notch  was  rectangular,  with  the  depth 
approximately  half  the  surface  length,  however,  the  crack  growing  from  the  notch  is  assumed  to 
quickly  attain  a  semi-elliptical  shape. 

A  Goodman  Diagram  approach  [37]  was  used  to  obtain  a  rough  estimate  of  the  equivalent 
maximum  applied  stress  levels  (amax)  at  various  R  values.  The  Goodman  line  or  the  constant 
lifetime  line  is  shown  in  Fig.  4  where  amax  of  860  MPa  at  R  =  0.05  is  taken  as  the  baseline 
condition.  It  is  recognized  that  the  Goodman  line  as  shown  here  is  only  an  approximation  of  a 
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constant-lifetime  line,  and  several  modifications  to  this  model  have  been  proposed  [38]. 
However,  for  the  purposes  of  determining  equivalent  stress  levels,  and  appropriate  combinations 
of  R,  applied  stress  level  and  notch  size  for  successful  small-crack  growth  tests,  this  approach 
can  be  used  to  obtain  reasonable  estimates.  The  intersection  points  between  the  Goodman  line 
and  lines  corresponding  to  R  =  -0.5  and  0.5  give  the  approximate  equivalent  applied  stress  levels 
at  these  stress  ratios.  The  corresponding  amax  values  applied  in  the  small-crack  growth  tests  are 
given  in  Table  I. 

Using  R  =  0.05  as  the  baseline,  a  notch  depth  of  5  pm  was  sufficient  to  produce  crack 
initiation  at  a  majority  of  notches  at  cjmax  =  860  MPa.  The  corresponding  initial  stress  intensity 
factor  range  (AK)  for  a  semi-elliptical  crack  of  the  same  depth  is  about  2.13  MPaVm,  as  given  in 
Table  I.  The  desired  initial  crack  depths,  needed  to  produce  the  initial  AK  of  2.13  MPaVm  at  R  = 
-0.5  and  0.5,  are  also  listed  in  Table  I.  Based  on  this  reasoning,  notch  depths  of  5  and  15  pm 
(with  the  surface  length,  2c,  being  twice  the  depth)  were  utilized  at  R  =  -0.5  and  0.5  respectively. 
As  indicated  in  Table  I,  small  cracks  initiated  from  a  majority  of  micro-notches  under  each  of  the 
three  combinations  of  R,  amax,  and  initial  notch  size. 

The  small-crack  growth  was  monitored  by  the  standard  acetate  replication  method  [39], 
Tests  were  periodically  interrupted  and  the  specimen  was  cooled  to  room  temperature  for 
replication.  Replicas  were  obtained  at  60%  of  the  maximum  load  in  a  fatigue  cycle.  An  optical 
microscope  was  used  to  measure  crack  lengths  from  replicas. 

Heating  system:  As  stated  previously,  the  small-crack  growth  tests  were  conducted  at 
260°C.  A  collet  gripping  mechanism  was  used  with  MTS  hydraulic  grips  to  load  the  specimens. 
A  four-zone  quartz-lamp  heating  system  was  employed  to  stabilize  specimens  at  the  test 
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temperature.  Temperatures  were  independently  controlled  and  monitored  at  four  points  along  the 
gauge  length  and  two  points  outside  the  gauge  on  both  ends  of  the  specimen.  Loop 
thermocouples,  that  were  spring  loaded  to  maintain  contact  with  the  specimen,  were  used  in  the 
gauge  section.  Figure  5  represents  the  thermocouple  configuration  in  a  dummy  specimen  which 
was  used  to  verify  the  heating  system.  In  the  dummy  specimen,  the  loop  thermocouples  were 
placed  in  the  same  configuration  as  in  an  actual  test,  but  extra  thermocouples  were  welded  to  the 
“side”  and  “front”  of  the  specimen,  as  illustrated  by  arrows  in  the  figure.  A  plot  of  the  stabilized 
temperature  profile,  i.e.,  after  the  ramp-up  period,  in  the  gauge  section  as  recorded  by  these 
thermocouples  is  shown  in  Fig.  6.  The  temperatures  were  recorded  for  duration  typical  of  a 
fatigue  test  and  found  to  be  very  stable  at  each  measurement  location  throughout  the  entire  time. 
As  shown  in  the  figure,  the  temperature  readings  were  within  ±3°C  of  the  set  point,  and  the  loop 
thermocouple  readings  were  in  close  agreement  with  those  of  the  welded  thermocouples. 


3.  SMALL-CRACK  GROWTH  DATA  REDUCTION 

3.1  K-Solution 

The  K-solution  by  Raju  and  Newman  [40]  for  a  semi-elliptical  surface  crack  in  a  rod  was 
used  in  this  study.  The  solution  at  the  surface  tip  of  the  crack  was  employed,  which  is  given  by 
[40]: 

K  =  S  t^-F  (1) 

yj  Q 

where  S  is  the  applied  stress  level,  a  is  the  half  crack  length,  Q  is  the  shape  factor  for  an  ellipse, 
and  F  is  the  boundary  correction  factor  [40],  The  crack  aspect  ratio,  a/c  of  1  was  assumed  for 
which  Q  =  2.464.  The  values  of  F  at  the  surface  tip  were  tabulated  in  [40]  as  a  function  of 
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normalized  crack  size,  a/D  where  D  is  the  rod  diameter.  A  2nd  order  polynomial  fit  to  this  tabular 
data  was  used  as  a  representation  of  the  factor,  F,  which  was:  F(j)  =  1.445  +  0.0727(j)  + 
2.5016.  While  for  positive  R,  AK  was  given  by  the  Maximum  K  -  Minimum  K  (Kmax  -  Kmin)  in 
a  fatigue  cycle,  for  negative  R,  only  the  positive  portion  of  the  loading  cycle  was  included  in  AK, 
i.e.,  AK  was  taken  as  Kmax. 

3.2  Calculation  of  growth  rates 

Typically,  the  small-crack  growth  behavior  shows  more  variability  than  the  long-crack 
regime,  due  to  stronger  influence  of  interactions  of  the  crack  front  with  the  local  microstructure 
[21,  22],  This  variability  is  reflected  both  in  terms  of  differences  in  growth  rates  between  cracks, 
as  well  as  fluctuations  in  growth  rate  of  a  given  crack.  Since  the  small  cracks  are  on  the  order  of 
the  microstructural  size  scale,  and  extension  of  a  crack  in  a  given  fatigue  loading  interval  can  be 
on  the  sub-micrometer  to  the  micrometer  scale,  the  accuracy  of  measurement  together  with  the 
frequency  of  data  collection  have  a  strong  bearing  on  the  calculated  growth  rate,  da/dN. 
Therefore,  in  order  to  isolate  the  microstructural  influence  on  the  small  crack  behavior,  it  is 
important  to  minimize  the  contributions  of  external  factors  on  the  fluctuation  in  growth  rate  in 
the  da/dN  vs.  AK  space.  Towards  this  end,  a  data  reduction  method  developed  by  Larsen,  et  al. 
[41]  was  used  in  this  study.  This  technique  will  be  referred  to  as  the  crack-length  increment 
method  in  the  remainder  of  this  paper.  The  method  is  based  on  a  polynomial  fit  to  crack  length 
(a)  vs.  cycles  (N)  data  blocks,  where  the  number  of  data  points  in  a  block  is  varied  to  satisfy  an 
error  criterion  [41],  The  data  block  is  moved  by  a  prescribed  crack  length  increment  (as  opposed 
to  point-wise  increment),  which  was  taken  to  be  on  the  order  of  the  microstructural  size  scale  (5 
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pm,  which  is  on  the  order  of  the  ap  size)  in  the  present  study.  An  illustration  of  data  reduced  by 
this  method  is  given  in  Fig.  7  where  the  crack-length  increment  method  is  compared  to  the  3- 
point  sliding  polynomial  method.  The  data  shown  in  the  figure  is  for  Ti-6-2-4-6  at  the 
temperature  of  260°C  and  R  =  -0.5.  Note  that  16  small  cracks,  which  initiated  from  FIB 
machined  micro-notches  in  two  specimens,  are  represented  in  Fig.  7. 


4.  RESULTS  AND  DISCUSSION 

4.1  Effect  of  Stress-Ratio  on  the  Small-Crack  Growth  Behavior 

The  effect  of  R  on  the  small  crack  growth  behavior  at  260°C  is  shown  in  Fig.  8  in  terms 
of  da/dN  vs.  the  nominal  AK.  As  shown,  R  of  -0.5,  0.05,  and  0.5  were  employed.  The  maximum 
stress  levels  applied  at  these  R  are  given  in  Table  1.  All  small  cracks  represented  in  the  figure 
initiated  from  FIB  machined  micro-notches.  Numbers  of  such  small  cracks  were  16  (in  two 
specimens),  7,  and  9  at  R  of  -0.5,  0.05,  and  0.5,  respectively  (Fig.  8).  A  qualitative  examination 
of  Fig.  8  suggests  that  R  did  not  produce  a  significant  influence  on  the  small-crack  growth 
behavior.  Similar  results  were  obtained  in  a  separate  study  on  the  effect  of  R  on  the  small-crack 
growth  regime  in  another  a+p  titanium  alloy,  Ti-6-4  [32],  Given  that  crack-growth  lifetime  can 
be  very  sensitive  to  the  statistical  parameters  of  the  small-crack  growth  behavior  [42],  from  the 
perspective  of  probabilistic  life  prediction,  the  next  questions  are:  (i)  how  to  probabilistically 
represent  the  small-crack  growth  data  for  a  life  prediction  analysis,  (ii)  what  is  the  statistical 
significance  of  the  differences  in  small-crack  behavior  with  respect  to  R,  and  (iii)  what  is  the 
impact  of  differences  in  the  statistical  parameters  on  the  probabilistic  life  prediction.  These  three 
factors  are  discussed  in  the  following  subsections. 
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4.2  Probabilistic  Representation  of  the  Small-Crack  Growth  Data 

Two  approaches  to  probabilistic  modeling  of  small-crack  growth  behavior  were  taken 
and  are  illustrated  here.  The  first  approach  is  based  on  fitting  a  model  to  the  reduced  da/dN  -  AK 
data,  which  requires  care  in  terms  of  verifying  if  the  model  accurately  produces  the  measured 
crack  length  vs.  cycles  ( a  vs.  N)  data.  The  second  approach  is  based  on  using  an  optimization 
method  to  calculate  the  model  parameters  that  give  the  best  match  between  the  predicted  and 
measured  a  vs.  N  data.  Both  methods  are  compared  in  the  following  subsections,  and  some 
considerations  that  were  found  to  be  important  in  incorporating  small-crack  growth  data  in  life 
prediction  are  highlighted.  To  study  the  effect  of  R  on  probabilistic  life  prediction,  the  latter 
method,  i.e.,  based  on  optimization  of  model  parameters,  was  implemented. 

Method  I  -  Fitting  a  model  to  the  reduced  da/dN  -  AK  data 

Choice  of  the  data  reduction  method  is  significant  if  the  resulting  da/dN  -  AK  data  is 
employed  in  fitting  a  model.  Particularly  in  probabilistic  life  prediction  where,  in  addition  to 
accurately  representing  individual  cracks,  the  variability  in  the  small-crack  growth  behavior  is 
quantitatively  represented,  minimizing  the  extraneous  contributions  to  the  da/dN  data  assumes 
added  significance.  In  the  method  of  fitting  a  model  to  the  reduced  da/dN-AK  data,  it  is  useful  to 
note  that  any  artifacts  in  the  a  vs.  N  measurements  may  be  further  enhanced  in  the  da/dN  -  AK 
space  as  discussed  in  a  recent  study  [43],  The  resulting  error  is  reflected  in  the  fit,  which,  due  to 
the  sensitivity  of  the  crack-growth  lifetime  to  the  da/dN  -  AK  representation  [42],  will  lead  to  an 
error  in  life  prediction.  The  crack-length  increment  method  by  Larsen,  et  al  [41],  described 
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previously,  was  employed  here,  which  seeks  to  address  the  contribution  of  artifacts  in  the  a  vs.  N 
data  in  the  reduced  da/dN  -  AK  data. 

One  objective  of  the  small-crack  growth  data  reduction  is  to  determine  the  model  that 
represents  the  da/dN  vs.  AK  relationship.  As  evident  from  Fig.  8,  a  power  law  is  justifiable  in 
this  case  as,  on  an  average,  the  small  cracks  show  a  linear  trend  in  the  Log(da/dN)  vs.  Log(AK) 
space.  Therefore,  a  Power-law  equation,  similar  in  form  to  Paris-law  [44],  was  fit  to  data  for 
individual  cracks  and  can  be  expressed  as: 


HI  =  eCM"  (2) 

A  set  of  values  for  C  and  n  were  thus  obtained  for  each  R.  The  distribution  in  C  and  n  can  be 
approximated  by  a  joint-normal  probability  density  where  these  two  random  variables  are 
strongly  negatively  correlated  [45].  Some  details  of  this  approach  as  applied  in  a  different  study 
are  given  in  [35],  The  method  is  illustrated  in  Fig.  9  with  the  help  of  the  set  of  values  for  C  and  n 
at  R  =  -0.5.  Figure  9(a)  shows  the  small-crack  growth  curves  reduced  by  the  crack-length  sliding 
method  at  R  =  -0.5  and  the  corresponding  power-law  fits.  Note  that  only  7  curves  are  shown  in 
Fig.  9(a)  although  16  of  the  18  FIB  notches  (in  two  specimens)  initiated  a  crack.  This  down- 
selection  pertains  to  a  final  crack  length  criterion  due  to  the  consideration  of  minimizing  the 
artificial  contributions  to  the  statistical  parameters  of  small-crack  growth  and  is  explained  later. 
The  normal-distribution  parameters  for  C  and  n,  based  on  these  7  curves,  were  derived  using  the 
maximum  likelihood  estimate.  The  joint-normal  density  representing  the  two  variables  is  given 
in  Fig.  9(b).  The  correlation  coefficient,  p  in  this  case  was  -0.97. 
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Modeling  individual  small-crack  growth  curves  by  a  power-law  is  expected  to  capture  the 
average  growth  behavior  of  a  given  crack  in  terms  of  da/dN  vs.  AK,  although  the  growth  rate 
fluctuates  about  the  average  due  to  microstructural  influences.  This  representation  may  be 
justified  if  the  calculated  a  vs.  N  based  on  the  modeled  behavior  reproduces  the  measured  a  vs. 
N  data.  Additionally,  in  order  to  verify  the  accuracy  of  the  present  approach,  it  is  important  to 
confirm  that  the  predicted  crack-growth  lifetime  between  a  given  crack-length  interval  matches 
the  experimentally  observed  change  in  lifetime  for  the  same  interval  as  suggested  by  McClung 
[46],  These  factors  are  demonstrated  in  Fig.  10.  The  calculated  a  vs.  N  based  on  a  power-law  fit 
to  the  da/dN  vs.  AK  data  using  the  crack-length  increment  method  of  data  reduction  for  a  crack  at 
R  =  -0.5  is  compared  to  the  measured  a  vs.  N  data  for  the  same  crack  in  Fig.  10(a).  A  reasonable 
agreement  between  the  calculated  and  measured  a  vs.  N  data  can  be  seen  from  this  figure.  Next, 
crack  growth  lifetimes  were  calculated  for  the  interval  between  the  first  measured  crack  length 
(i.e.,  when  a  crack  was  first  detected  to  have  grown  out  of  the  FIB  notches)  and  the  final  length 
for  the  given  crack.  The  calculated  vs.  the  measured  crack-growth  lifetime  for  all  cracks  included 
in  the  probabilistic  model  at  R  =  -0.5  are  plotted  in  Fig.  10(b).  As  shown,  the  predicted  crack 
growth  lifetimes  were  within  10%  of  the  measured  lifetimes  between  the  respective  crack  growth 
intervals.  This  suggests  that  this  method  represents  a  reasonable  treatment  of  small-crack  data  for 
the  purpose  of  probabilistic  life  prediction. 

Considerations  in  fitting  a  model  to  the  da/dN  —  AK  data :  Variability  in  the  small-crack 
growth  rates,  both  in  terms  of  the  fluctuations  in  growth  rate  of  a  given  crack  and  differences 
between  cracks  growing  in  nominally  the  same  material,  is  expected  due  to  a  strong  interaction 
of  small  cracks  with  the  local  micro  structure  and  the  fact  that  subsurface  crack  growth  may 
affect  the  path  and  velocity  of  the  crack  as  measured  on  the  specimen  surface.  The  goal  here  is  to 
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minimize  artificial  contributions  to  this  variability.  As  alluded  to  previously,  error  in  the  da/dN  - 
AK  data  can  arise  from  external  factors  such  as  data  collection  frequency  and  measurement 
resolution  or  measurement  error.  In  probabilistic  representation,  an  additional  factor  that  was 
found  important  was  the  final  crack  length.  The  first  two  factors  listed  above  are  interrelated,  in 
that  the  influence  of  measurement  resolution  or  error  will  get  exacerbated  with  higher  data 
collection  frequency,  and  is  attributed  to  the  fact  that  the  increment  in  the  length  of  a 
microstructurally-small  crack  can  be  on  the  order  of  the  resolution  of  optical  measurement  or 
within  the  range  of  experimental  error.  These  variables  artificially  contribute  to  the  fluctuations 
in  the  reduced  da/dN  vs.  AK  data.  The  third  factor  listed  above,  i.e.,  the  final  crack  length  does 
not  affect  the  fluctuation  in  growth  rates  of  the  given  crack  but  may  artificially  distort  the 
statistical  parameters  of  a  small-crack  growth  model  if  the  model  is  to  be  used  in  prediction  of 
the  total  crack  growth  lifetime.  These  factors  are  illustrated  in  a  following  discussion. 

An  illustration  of  the  effect  of  the  external  contributions,  listed  previously,  to  the  reduced 
da/dN  -  AK  data  on  the  probabilistic  representation  of  small-crack  growth  is  given  in  Fig.  1 1  (a) 
and  (b)  with  the  help  of  the  small  cracks  measured  at  R  =  -0.5.  The  data  collection  interval  at  the 
beginning  of  this  test  was  500  cycles,  which  was  increased  to  1000  cycles  in  the  intermediate 
growth-rate  period,  and  then  decreased  back  to  500  cycles  towards  the  end  of  the  test.  In  Fig. 
11(a),  data  reduced  by  the  crack-length  increment  method  are  compared  to  the  3-point  and  5- 
point  sliding  polynomial  methods.  Power-law  representations  of  the  da/dN  -  AK  data  are  also 
shown  for  each  method.  In  case  of  the  sliding  polynomial  method,  the  variables,  C  and  n  depend 
on  the  frequency  of  small-crack  measurements  and  the  number  of  points  in  the  polynomial  fit.  In 
this  example  (Fig.  11(a)),  a  fit  to  the  3-point  sliding  polynomial  data  misrepresents  the  average 
trend  in  da/dN  with  respect  to  AK  and  gives  values  of  C  and  n  that  are  significantly  different  than 
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those  given  by  the  power-law  fit  of  data  by  the  crack-length  increment  method.  A  power-law  fit 
to  the  data  by  the  5-point  sliding  polynomial  method  gives  a  result  closer  to  that  by  the  crack- 
length  increment  method.  The  impact  of  error  in  the  power-law  variables,  C  and  n  can  be 
illustrated  by  comparing  the  measured  vs.  predicted  a  vs.  N  data,  which  is  presented  in  Fig.  11 
(b)  .  In  this  case,  the  3-point  sliding  method  produces  significant  errors  in  the  predicted  a  vs.  N 
which  will  be  reflected  in  the  calculated  crack-growth  lifetimes.  On  the  other  hand,  the  5-point 
sliding  polynomial  compares  well  with  the  crack-length  increment  method. 

As  mentioned  before,  the  above  factors  may  be  further  accentuated  in  a  probabilistic 
model  of  the  small-crack  growth  behavior  when  a  set  of  cracks  at  a  given  condition  is 
considered.  To  illustrate  this  point,  power-law  fits  to  the  small-crack  data  set  at  R  =  -0.5,  where 
data  were  reduced  by  the  3-point  sliding  polynomial  method  are  compared  to  power-law 
representations  of  the  da/dN  -  AK  data  by  the  crack-length  increment  method  in  Fig.  1 1(c).  The 
statistical  parameters  of  the  model  variables,  listed  in  the  figure,  demonstrate  the  strong  influence 
of  the  method  on  these  parameters.  For  example,  the  mean  and  standard  deviation  of  C  by  the  3- 
point  sliding  polynomial  method  were  -24.33  and  3.55  respectively.  On  the  other  hand,  those 
parameters  were  -22.10  and  0.74  when  the  crack-length  sliding  method  was  used.  This  difference 
in  statistical  parameters  of  C  will  likely  have  a  significant  effect  on  the  predicted  lifetime 
distribution,  as  evident  from  Fig.  11  (b). 

The  statistical  parameters  of  the  small-crack  growth  variables  in  a  probabilistic 
representation  were  also  found  to  be  sensitive  to  the  final  crack  length  measurement  for  a  given 
crack.  This  is  because,  if  a  crack  did  not  grow  to  a  sufficient  size,  the  power-law  fit  to  the 
corresponding  da/dN  -  AK  data  may  be  biased  by  a  local  snapshot  of  the  growth  behavior,  which 
can  be  very  different  from  the  overall  trend  if  that  crack  had  grown  to  a  larger  size.  This  factor 
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can  artificially  distort  the  statistical  parameters  of  the  small-crack  growth  model  variables  if  the 
model  is  to  be  used  in  calculation  of  the  distribution  in  the  total  crack-growth  lifetime,  i.e.,  until 
fracture.  This  is  illustrated  in  Fig.  12,  where  the  same  crack  as  in  Fig.  11(a)  is  shown.  The  final 
measurement  for  this  crack  was  2c  =  928  pm,  which  corresponded  to  a  AK  of  about  23  MPa- 
m1/2.  The  power-law  fit  to  the  da/dN  vs.  AK  data  is  shown  in  the  figure.  For  the  sake  of 
argument,  if  other  cracks  in  the  specimen  were  more  dominant,  then  the  final  measurement  for 
this  crack  just  before  fracture  would  have  been  smaller.  To  illustrate  the  point,  the  da/dN  -  AK 
data  assuming  that  the  crack  length  at  the  time  of  failure  was  only  70  pm  (AK  ~  6  MPa-m  )  is 
shown  in  Fig.  12.  The  corresponding  power-law  fit  is  very  different  from  when  the  crack  grew  to 
a  larger  length,  as  it  is  skewed  by  a  local  trend,  which  is  not  representative  of  the  true  behavior 
across  the  full  range  of  crack  length.  The  values  of  C  for  the  full  and  the  truncated  crack  were  - 
23.39  and  -21.80  respectively  (Fig.  12),  which  is  a  significant  difference  in  terms  of  lifetime 
calculation.  The  error  due  to  this  factor  can  be  minimized  if  cracks  smaller  than  a  critical  size  are 
not  included  in  the  probabilistic  representation  of  the  small-crack  growth  regime.  Ideally,  the 
critical  crack  length  should  correspond  to  the  AK  beyond  which  a  fit  to  the  data  is  largely 
independent  of  the  final  crack  length.  The  crack  size  corresponding  to  the  transition  to  a  large 
crack  may  be  a  reasonable  critical  length,  as  it  can  be  suggested  that  once  a  crack  has 
transitioned  to  a  large  crack,  its  overall  trend  is  already  developed  and  a  model  representation  of 
the  data  will  likely  be  minimally  influenced  by  further  increase  in  length  beyond  the  transition 
point.  This  criterion  was  adopted  here,  and  only  cracks  that  had  grown  to  the  half  crack  length  of 
at  least  100  pm  (which  is  about  20  times  the  ap  size)  at  the  time  of  specimen  fracture  were 
considered  in  the  probabilistic  representation.  This  is  based  on  the  study  by  Lankford,  et  al  who 
reported  that  in  many  materials  the  transition  from  the  small  to  large-crack  behavior 
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corresponded  to  about  10-20  times  the  controlling  microstructural  unit  size  [47].  It  should  be 
noted  that,  while  the  final  crack  length  criterion  is  necessary  to  obtain  accurate  statistical 
parameters  in  calculation  of  total  small-crack  growth  lifetimes,  this  may  bias  the  population  of 
cracks  used  in  the  analysis  towards  those  that  exhibited  faster  growth.  This  is  an  artifact  of 
multiple  propagating  cracks  in  each  specimen  due  to  multiple  FIB  notches,  which  suppresses  the 
full  a  vs.  N  information  from  some  slow  growing  cracks.  An  ideal  case  will  be  when  each 
specimen  has  single  notch  that  initiates  the  only  propagating  crack.  Then,  if  a  sufficient  number 
of  tests  are  conducted,  a  more  complete  range  of  growth  rates  can  be  represented  in  the  model  of 
small-crack  growth.  In  terms  of  probabilistic  life  prediction,  the  influence  of  this  artifact  will  be 
to  produce  predictions  that  are  conservative  with  respect  to  the  actual  distributions. 

It  is  useful  to  note  that,  in  many  cases,  a  standard  data  reduction  method  such  as  sliding 
polynomial  may  be  sufficient  in  probabilistic  life  prediction  but  the  reduced  data  will  be 
dependent  upon  the  frequency  of  measurements  and  the  number  of  points  in  a  polynomial  fit, 
which  will  influence  the  statistical  parameters  in  a  probabilistic  representation.  In  some  cases 
(depending  on  the  amount  of  data),  a  higher-point  sliding  polynomial  method  may  yield  a  similar 
result  as  the  crack-length  increment  method  used  here,  as  demonstrated  by  the  above  example. 
However,  an  important  feature  of  the  latter  method  is  that  it  provides  a  consistent  da/dN  -  AK 
data  independent  of  experimental  variables  such  as  data  collection  frequency  and  measurement 
error  as  long  as  an  appropriate  crack-length  increment  is  adopted.  As  discussed  above,  a  crucial 
metric  for  the  accuracy  of  a  small-crack  growth  model  is  to  compare  the  a  vs.  N  behavior 
predicted  by  the  model  representation  to  the  measured  a  vs.  N  data  [46].  It  can  be  said  that,  to  the 
extent  that  a  model  limits  the  error  in  the  predicted  a  vs.  N  to  an  acceptable  level,  it  might  be 
applicable  in  probabilistic  life  prediction  analysis. 
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Method  II  -  Method  of  optimization  of  the  small-crack  growth  model  parameters 

As  illustrated  above,  fitting  a  model  to  the  da/dN  -  AK  data  by  the  crack-length 
increment  method  produced  an  acceptable  level  of  error  in  calculated  crack-growth  lifetimes 
(Fig.  10).  However,  a  more  rigorous  approach  can  be  to  calculate  the  optimal  model  parameters 
that  will  give  the  best  match  between  the  predicted  and  the  measured  a  vs.  N  data.  This  method 
can  further  improve  the  agreement  between  the  predicted  and  actual  crack  growth  lifetimes  since 
the  objective  is  to  accurately  represent  the  measured  a  vs.  N  data  rather  than  fitting  to  the 
reduced  da/dN  -  AK  data,  which  is  sensitive  to  measurement  errors. 

The  data  reduced  by  the  crack  length  increment  method  indicated  that  a  power-law  model 
(Eqn.  2)  was  a  reasonable  representation  of  the  small-crack  growth  behavior.  The  optimization 
problem,  in  this  case,  is  to  find  the  optimal  values  of  the  power-law  variables,  C  and  n  that  will 
minimize  the  error  between  the  predicted  and  the  measured  a  vs.  N  data.  The  cost  function,  f(C, 
n)  that  has  to  be  minimized  was  formulated  as  the  least  squares  error  between  the  predicted  and 
the  measured  a  vs.  N  data,  given  by: 

f{C,Tl)  —  Yii-N^^precL  ~  ameas 0) 

Where,  aPred  and  ameas  are  the  predicted  and  the  measured  half  crack  length,  respectively.  Ni 
corresponds  to  the  cycle  count  at  which  the  crack  was  first  detected  and  Nn  represents  the  cycle 
count  at  the  final  measurement.  The  summation  is  performed  over  all  cycle  counts,  starting  with 
Ni,  at  which  crack  length  measurements  were  recorded.  Note  that  a  weighted  least  squares  error 
[48]  can  also  be  explored  as  a  formulation  for  the  cost  function,  which  may  be  more  appropriate 
when  the  data  collection  frequency  is  high  and  the  data  exhibits  non-uniform  variance  with 
respect  to  the  cycle  count.  A  Mathematica™  script  was  written  to  perform  the  optimization 
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calculation  in  which  the  Nelder-Mead  numerical  optimization  method  was  implemented.  The 
range  of  values  of  C  and  n  obtained  by  the  first  approach,  i.e.,  fitting  the  da/dN  -  AK  data,  were 
used  as  the  initial  variable  search  space  in  the  optimization  routine.  The  a  vs.  N  predictions  for  a 
given  C  and  n  were  made  by  the  cycle-by-cycle  calculation  of  crack  extension. 

An  example  of  the  predicted  a  vs.  N  data  using  the  values  of  C  and  n  given  by  the 
optimization  method  is  shown  in  Fig.  13(a)  for  a  crack  at  R  =  -0.5  (which  is  the  same  crack  as 
shown  in  Fig.  10(a)).  For  comparison,  the  predicted  a  vs.  N  using  the  C  and  n  values  from  fitting 
a  power-law  to  the  da/dN  -  AK  data  is  also  included.  It  is  not  surprising  that  the  optimized  model 
produces  a  better  agreement  between  the  predicted  and  the  measured  a  vs.  N  data.  The  power- 
law  representations  of  small-crack  curves  at  R  =  -0.5  based  on  the  optimal  C  and  n  values  are 
compared  to  those  based  on  the  first  approach  in  Fig.  13(b).  As  shown,  the  model  based  on 
optimized  variables  does  not  differ  significantly  from  the  power-law  fits  to  data  reduced  by  the 
crack-length  increment  method,  indicating  that  only  small  changes  in  the  values  given  by  the 
first  approach  are  required  to  attain  optimum  variables.  In  Fig.  13(c),  the  actual  vs.  predicted 
small-crack  growth  lifetimes  using  the  optimized  C  and  n  are  plotted.  As  expected,  an  improved 
agreement  between  the  actual  and  predicted  values  is  achieved  when  compared  to  Fig.  10(b).  As 
shown,  in  this  case,  the  predictions  by  the  optimization  method  were  within  ±2%  of  the  actual 
values  compared  to  ±10%  by  the  first  method.  As  stated  before,  the  distributions  in  C  and  n  at 
various  R  as  given  by  the  optimization  method  were  used  in  the  subsequent  probabilistic  life 
prediction  analysis. 

4.3  Statistical  Significance 
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To  determine  the  statistical  significance  of  the  effect  of  R  on  the  small-crack  growth 
behavior,  a  t-test  [49]  was  employed,  which  is  a  standard  technique  to  assess  if  two  data  sets 
have  statistically  different  means.  Since  the  mean  and  variance  of  the  simulated  crack-growth 
lifetime  distribution  were  found  to  be  most  sensitive  to  C  [42],  the  test  was  performed  on  this 
random  variable.  Only  the  populations  of  cracks  at  each  R  that  were  employed  in  the 
probabilistic  representation  of  small  crack  growth  were  tested.  A  two  sample  t-test  assuming 
unequal  variances  [49]  was  used.  The  results  of  this  test  are  given  in  Table  II.  The  t  value  for  the 
null  hypothesis  between  data  pairs  at  R  =  0.05  and  0.5,  R  =  0.05  and  -0.5,  and  R  =  0.5  and  -0.5 
were  -2.19,  -2.41,  and  -0.59  respectively  (Table  II).  The  corresponding  two-tail  P  values  (i.e.,  the 
probability  that  the  t  value  is  less  than  the  critical  t)  were  0.07,  0.04,  and  0.56  respectively.  One 
interpretation  of  this  test  can  be,  that  there  is  low  probability  that  the  populations  of  C  at  R  = 
0.05  and  0.5  and  those  at  R  =  0.05  and  -0.5  have  statistically  the  same  mean.  Also,  there  is 
moderate  probability  of  the  populations  at  R  =  0.5  and  -0.5  to  have  statistically  different  means. 
Given  that  the  calculated  lifetime  distribution  is  very  sensitive  to  the  mean  and  standard 
deviation  of  C,  the  next  question  is  whether  these  differences  are  significant  in  terms  of  life 
prediction,  as  is  discussed  in  the  next  sub-section. 

4.4  Effect  on  Probabilistic  Life  Prediction 

The  Life  Prediction  Method 

Previous  studies  on  life-limiting  fatigue  behavior  in  several  structural  materials  have 
shown  that  the  lifetime  distribution  separates  into  two  populations,  termed  as  mean-dominating 
and  life-limiting  populations  [29,  35,  36],  The  life-limiting  population  is  defined  as  the  short¬ 
lifetime  population  that  follows  a  different  trend  than  the  mean-dominating  population  with 
respect  to  stress  level  or  any  other  variable.  More  specifically,  the  mean  behavior  diverged  (or 
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converged)  with  respect  to  the  life-limiting  behavior  as  a  function  of  microstructure,  stress  level, 
temperature,  etc.  [29,  35,  36].  For  the  probability  of  occurrence  of  the  life-limiting  failures  varied 
as  a  function  of  these  variables,  for  example,  this  probability  decreased  with  a  decrease  in  the 
stress  level  [36].  Through  probabilistic  simulation  of  the  crack  growth  lifetime,  it  was  shown  that 
the  distribution  in  the  life-limiting  points  can  be  modeled  by  the  small  +  long  crack  growth 
lifetime  distribution  for  cracks  initiating  from  the  controlling  microstructural  scale  [35,  36].  The 
simulation  was  performed  by  the  Monte  Carlo  method  where  the  input  random  variables  were: 
(i)  the  crack-initiation  size,  (ii)  the  small-crack  growth  variable,  C,  and  (iii)  and  the  small-crack 
growth  variable,  n. 

The  crack-initiation  size  was  modeled  by  the  lognormal  distribution,  the  parameters  of 
which  were  derived  by  maximum  likelihood  estimate  (MLE)  based  on  a  population  of  crack- 
initiation  sizes  measured  in  fractured  fatigue  specimens.  For  example,  a  typical  crack  initiation 
site  in  duplex  Ti-6-2-4-6  at  260°C  is  shown  in  Fig.  14  (a)  [50].  Crack  initiation  in  this  alloy  is 
accomplished  by  formation  of  a  crystallographic  facet  across  a  (or  a  few)  primary-a  particle,  as 
indicated  in  the  figure.  While  subsurface  failures  were  observed  under  certain  loading  conditions, 
the  life-limiting  failures  were  always  associated  with  surface  crack  initiation.  Areas  of  the  crack- 
initiation  facets  were  measured  in  a  number  of  fractured  specimens  and  the  distribution  in  the 
equivalent  crack-initiation  depth,  a  (i.e.,  radius  of  a  semicircle  with  the  same  area  as  the  crack- 
initiation  facet(s))  at  R  =  0.05  is  shown  in  Fig.  14(b)  in  lognormal  cumulative  distribution 
function  (CDF)  space.  The  corresponding  lognormal  fit  to  the  data,  which  was  used  as  an  input 
to  the  probabilistic  model,  is  also  shown.  The  crack-initiation  size  was  assumed  to  not  vary  with 
respect  to  R,  which  will  be  the  case  if,  in  life-limiting  specimens,  the  crack-initiation  mechanism 
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involved  facet  formation  in  primary-a  particles  at  the  surface.  However,  this  was  not  verified  at 
all  stress  ratios. 

The  two  random  variables  associated  with  the  small-crack  growth,  C  and  n,  were 
described  previously.  A  set  of  values  of  C  and  n  were  obtained  at  each  R,  which  corresponded  to 
the  optimal  values  given  by  the  method  of  optimization  of  the  small-crack  growth  model  with 
respect  to  the  measured  a  vs.  N  data.  The  distribution  in  C  and  n  were  modeled  by  the  joint- 
normal  probability  distribution  function,  the  parameters  of  which  were  derived  my  maximum 
likelihood  estimate,  as  discussed  before.  The  joint  probability  densities  representing  C  and  n  at  R 
=  -0.5,  0.05,  and  0.5  that  were  used  in  the  probabilistic  model  are  shown  in  Fig.  15  and  the 
corresponding  parameters  are  given  in  Table  III.  Note  that  the  correlation  coefficient,  p,  which, 
in  this  case,  borne  from  the  optimized  representations  of  small-crack  growth  were  -0.97,  -0.99, 
and  -0.99  respectively  at  R  =  -0.5,  0.05,  and  0.5. 

From  the  perspective  of  predicting  a  limiting  fatigue  lifetime  for  design  purposes,  the 
above  model  provides  a  mechanism-based  method  to  account  for  the  worst-case  fatigue  failure 
distribution,  almost  independent  of  both  the  amount  of  fatigue  data  and  whether  a  life-limiting 
point  is  captured  by  the  lifetime  data  from  S-N  tests.  For  illustration,  in  Fig.  16,  the  predicted 
life-limiting  distribution  is  compared  to  the  experimental  data  for  the  case  of  R  =  0.05.  The 
predicted  distribution  consists  of  10,000  Monte  Carlo  predictions.  The  plot  represents  the  CDF 
based  on  the  lognormal  distribution.  As  shown,  the  experimental  points  display  a  step-like 
morphology  in  the  CDF  space,  indicating  a  superposition  of  two  distributions.  The  three  shortest 
lifetime  points,  which  form  the  lower  part  of  the  step,  are  considered  as  belonging  to  the  life- 
limiting  population.  The  predicted  distribution  had  a  similar  slope  to  the  life-limiting  points 
when  those  are  plotted  as  a  separate  distribution  (Fig.  16).  However,  although  the  predicted 
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distribution  covers  the  range  of  lifetimes  of  the  life-limiting  points,  it  is  conservative  with  respect 
to  the  average  behavior  of  the  life-limiting  population.  This  discrepancy  could  be  due  to  the 
factors  that  only  a  limited  number  of  life-limiting  observations  are  being  compared  to  the 
simulation  and  that  a  very  small  number  of  cracks  (4  cracks  that  grew  to  a  >  100  pm  at  R  =  0.05) 
were  used  to  define  the  input  distributions  of  the  small-crack  growth  random  variables.  Also,  as 
discussed  previously,  the  method  of  using  single  specimen  to  characterize  multiple  cracks  may 
introduce  a  bias  in  the  analysis  towards  faster  cracks  since  the  full  crack-growth  trend  from 
slower  cracks  in  the  same  specimen  are  suppressed.  Nevertheless,  the  model  appears  to  produce 
a  conservative  estimate,  within  1.5X  of  the  mean  of  the  life-limiting  distribution.  If  the 
probabilistic  lifetime  limit  is  taken  as  the  1  in  1000  (termed  as  the  B0.1  lifetime)  probability  of 
failure,  the  predicted  lifetime  limit  at  R  =  0.05  was  about  12,856  cycles,  as  indicated  in  the 
figure.  The  present  study  employed  this  probabilistic  life-prediction  method  with  a  goal  to 
quantify  the  importance  of  the  differences  in  the  small-crack  growth  behavior  as  a  function  of  R 
on  the  probabilistic  lifetime  limit. 

Effect  of  R  on  probabilistic  lifetime  limit 

The  method  described  above  was  used  to  predict  the  life-limiting  distributions  at  the  three 
R  values,  which  are  compared  in  Fig.  17.  As  stated  before,  the  same  crack-initiation  size 
distribution  was  used  for  all  Rs.  Since  the  predictions  also  include  the  effect  of  the  applied  stress 
amplitude  on  lifetime,  in  order  to  represent  the  true  effect  of  the  small-crack  growth  differences 
on  life  prediction,  three  cases  were  considered.  In  each  case,  a  baseline  loading  condition  was 
chosen,  while  the  effect  of  employing  the  small-crack  behavior  at  various  Rs  on  the  predicted 
distributions  was  studied.  In  other  words,  if  the  differences  between  the  small-crack  growth 
behaviors  with  respect  to  R  were  insignificant,  the  predicted  distributions  at  a  given  loading 
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condition  should  be  very  similar.  Figures  17  (a),  (b)  and  (c)  show  predictions  for  the  loading 
parameters  corresponding  to  the  applied  R  of  -0.5,  0.05,  and  0.5,  respectively.  As  shown,  in  each 
of  the  three  cases,  the  small-crack  data  employed  in  the  model  has  a  significant  impact  on  the 
predicted  distributions.  This  can  be  further  illustrated  by  the  B0.1  lifetimes,  which  are  indicated 
in  Figs.  17  (a)  to  (c),  and  show  differences  ranging  from  26%  to  72%  with  respect  to  the  baseline 
condition.  Therefore,  in  the  present  problem,  the  differences  in  the  statistical  parameters  of  the 
small-crack  growth  with  respect  to  R  are  significant  from  a  life-prediction  perspective. 

The  above  analysis  demonstrates  that,  although  a  subjective  examination  of  the  present 
small-crack  growth  data  suggests  an  insignificant  effect  of  R,  the  effect  cannot  be  ignored  in 
probabilistic  life  prediction  and  particularly  in  probabilistic  life  limits.  This  also  emphasizes  the 
need  for  ensuring  that  the  modeled  da/dN  -  AK  relationship  accurately  reflects  the  measured  a 
vs.  N  data.  It  can  be  said  that  the  small-crack  growth  is  a  very  important  domain  of  fatigue,  both 
from  the  consideration  of  predicting  life-limiting  failures  and  due  to  the  sensitivity  of  the 
predictions  to  the  small-crack  growth  parameters.  Therefore,  the  small-crack  regime  needs  to  be 
accounted  for  in  probabilistic  life  prediction  but  care  in  terms  of  model  representation  of  small- 
crack  growth  is  required.  In  this  regard,  the  method  of  optimizing  the  small-crack  growth  model 
with  respect  to  the  measured  a  vs.  N  data  may  provide  a  means  to  enable  the  incorporation  of  the 
small-crack  regime  in  life  prediction  models. 

5.  CONCLUSIONS 

The  following  main  conclusions  can  be  drawn  from  this  study: 

(1)  Small-crack  growth  is  a  crucial  regime  of  fatigue  in  probabilistic  life  prediction  and  plays 

a  key  role  in  the  distribution  in  lifetimes  of  the  life-limiting  population. 
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(2)  When  employing  small-crack  growth  data  in  probabilistic  life  prediction,  extraneous 
factors  that  may  introduce  artificial  contributions  to  the  statistical  parameters  of  the 
small-crack  growth  model  variables  should  be  minimized  and  the  model  representation  of 
the  da/dN  vs.  AK  relationship  must  be  evaluated  against  predicted  vs.  measured  a  vs.  N 
data. 

(3)  The  probabilistic  predictions  revealed  that,  although  the  effect  of  R  on  the  small-crack 
growth  behavior  appeared  qualitatively  insignificant,  the  impact  of  the  statistical 
differences  with  respect  to  R  was  significant  on  the  probabilistic  lifetime  limit.  The  BO.  1 
lifetime  varied  in  the  range  of  26  and  72%  when  the  effect  of  R  on  the  small-crack 
growth  was  included. 
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Tables: 


Table  I:  Estimation  of  FIB  notch  depths  and  maximum  stress  levels  at  various  R 


R 

^max 

(Goodman) 

(MPa) 

AKforthe 
notch  depth  of 

5  pm  (MPaVm) 

Notch  depth 
forAK=  2.13 
(pm) 

Notch  depth 
employed 
(pm) 

Fraction  of 
notches  that 
initiated  a  crack 

-0.5 

805 

2.11 

5.1 

5 

0.9 

0.05 

860 

2.13 

5.0 

5 

0.8 

0.5 

925 

1.20 

15.6 

15 

1 
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Table  II:  Results  of  the  Student’s  t-test  on  data  sets  of  C  for  various  R  pairs 


R  pairs 

Hypothesis 

t  stat 

P-value 

(two-tail) 

^critical 

0.05  and 
0.5 

Null 

-2.19 

0.07 

2.45 

0.05  and  - 
0.5 

Null 

-2.41 

0.04 

2.26 

0.5  and  - 
0.5 

Null 

-0.59 

0.56 

2.26 

Table  III:  Statistical  parameters  of  small-crack  growth  model  variables 
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Random 

variable 

Distribution 

R  =  -0.5 

R  =  0.05 

R  =  0.5 

ft 

a 

ft 

a 

a 

a 

Lognormal 

2.29 

0.27 

2.29 

0.27 

2.29 

0.27 

C 

Joint-normal* 

-22.06 

0.77 

-23.04 

0.44 

-22.29 

0.39 

n 

Joint-normal* 

2.08 

0.30 

2.31 

0.28 

2.11 

0.24 

*The  correlation  coefficient,  p  between  C  and  n  were  -0.97,  -0.99,  and  -0.99 
atR  =  -0.5,  0.05,  and  0.5  respectively 
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Figures: 


20pm 

Fig.  1:  A  back-scattered  electron  image  of  the  microstructure  of  the  Ti-6-2-4-6  alloy  employed 
this  study 
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Strain  (mm/mm) 

Fig.  2:  Tensile  behavior  of  Ti-6-2-4-6  at  260  °C 
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1  mm 


10  (j.m 


Fig.  3:  FIB  machined  micro-notch  at  a  fatigue  crack  origin  of  a  small-crack  growth  specimen 
tested  at  R  =  0.5;  (a)  an  overview  of  the  crack  that  initiated  from  the  FIB  micro-notch  at  the  time 
of  final  fracture  and  (b)  crack  origin. 
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a  (MPa) 

m 

Fig.  4:  Goodman  diagram  approach  to  estimate  the  maximum  stress  levels  at  various  R 


36 

Approved  for  public  release;  distribution  unlimited. 


Front  1 

Front 2 

Front  3 


Side  1 
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Fig.  5:  Schematic  illustration  of  the  locations  of  loop  thermocouples  (which  were  used  to  control 
the  temperature  in  an  actual  test)  and  welded  thermocouples  in  a  dummy  specimen  that  was  used 
to  test  the  temperature  profde 
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Fig.  6:  Temperature  readings  from  the  loop  and  welded  thermocouples  in  a  dummy  specimen 


38 

Approved  for  public  release;  distribution  unlimited. 


Fig.  7:  A  comparison  of  small  fatigue  crack  data  reduction  by  the  3-point  sliding  polynomial  and 
crack-length  increment  methods. 
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Fig.  8:  Effect  of  R  on  the  small  crack  growth  behavior  in  Ti-6-2-4-6  at  260°C 
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(b) 

Fig.  9:  Probabilistic  representation  of  the  small  fatigue  crack  growth  data;  (a)  Power-law  fits  to 
the  small-crack  growth  curves  and  (b)  Joint-normal  probability  density  representing  the 
distribution  in  the  random  variables  C  and  n. 
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(b) 

Fig.  10:  Testing  the  accuracy  of  small-crack  growth  representation  by  power-law  fit  to  data 
reduced  by  the  crack-length  increment  method;  (a)  predicted  vs.  measured  a  vs.  N  data  and  (b) 
actual  vs.  predicted  small-crack  growth  lifetimes. 
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Fig.  11: 
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Fig.  11:  Effect  of  the  small-crack  growth  data  reduction  method  on  the  statistical  parameters  of  C 
and  n;  (a)  Illustration  of  the  effect  of  the  method  on  the  power-law  fit,  (b)  Illustration  of  the  error 
in  the  predicted  a  vs.  N  data,  and  (c)  the  effect  on  the  statistical  parameters,  mean  (pi)  and 
standard  deviation  (a),  of  C  and  n. 
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Fig.  12:  Illustration  of  the  effect  of  the  crack  length  at  the  time  of  fracture  on  the  fit  to  the  da/dN 
vs.  AK  data. 
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Fig.  13: 
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(C) 

Fig.  13:  Method  of  optimization  of  the  small-crack  growth  model  with  respect  to  the  measured  a 
vs.  N  data;  (a)  Comparisons  of  measured  and  predicted  a  vs.  N  using  the  crack-length  increment 
data  reduction  approach  and  the  optimization  method,  (b)  Optimized  power-law  representations 
of  small-crack  growth  curves  compared  to  fits  to  da/dN  -  AK  data  reduced  using  the  crack-length 
increment  method,  and  (c)  actual  vs.  predicted  crack  growth  lifetimes  for  cracks  at  R  =  -0.5  using 
the  optimized  representations  of  small-crack  growth. 
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Fig.  14:  Crack-initiation  size  distribution  in  Ti-6-2-4-6  at  260°C;  (a)  a  crack-initiation  site 
showing  a  faceted  ap  particle  at  the  crack  origin  (amax  =  860  MPa,  R  =  0.05)  and  (b)  The  crack 
initiation  size  distribution  shown  in  the  lognormal  cumulative  distribution  function  space. 
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Fig.  15:  Joint  normal  distributions  of  C  and  n  at  R  =  -0.5,  0.05,  and  0.5  used  in  probabilistic  life 
prediction. 
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Fig.  16:  Predicted  life-limiting  failure  distribution  at  R  =  0.05  and  260°C. 
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Fig.  17: 
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Fig.  17:  Influence  of  the  effect  of  R  on  the  small-crack  growth  regime  on  probabilistic  life 
prediction.  Predicted  crack-growth  lifetime  distributions  and  B0.1  lifetimes  at  (a)  R  =  -0.5,  (b)  R 
=  0.05,  and  (c)  R  =  -0.5. 
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